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£NJ , Abstract. Extending the work of Jia and Sverak on self-similar solutions of the Navier- 

» 1 Stokes equations, we show the existence of large, forward, discrete self-similar solutions. 

O: 

O ! 1 Introduction 

Denote = K 3 x (0,oo). Consider the 3D incompressible Navier-Stokes equations for 
1 velocity u : Ri — > R 3 and pressure p : Ri — > R, 

< 

d t u - Au + (u • V)u + Vp = 0, divu = 0, (1.1) 

-i— > ■ 

g \ in M^_, coupled with the initial condition 

u\ t=0 = u , divn = 0. (1.2) 

7— I ■ 

The system (jl.ip enjoys a scaling property: If u(x, t) is a solution, then so is 
|>. u^(x,t) := Xu{Xx,X 2 t) (1.3) 

<n: 

• for any A > 0. We say u(x, t) is self-similar (SS) if u = for every A > 0. In that case, 

the value of u(x, t) is decided by its value at any time moment £ = 3- and 

1 

>: u(x,t) = X(t)U(X(t)x), \(t) = -=, (1.4) 

.!— 1. V lot 

X: 

where f/(x) = u(x, ^) and a > is a parameter. On the other hand, if u = only for 
one particular A > 1, we say u is discrete self-similar (DSS) with factor A, or A-DSS. Its 
value in R\ is decided by its value in the strip i£l 3 and 1 < t < A 2 . We consider being 
SS as a special case of being DSS, and would say "strictly DSS" to exclude the former. We 
call them forward to indicate they are defined for t > 0. We can also consider (jl.ip for 

1. (x,t) £l 3 x (-oo,0), or 

2. x £ ]R 3 , u = u(x) is time independent. 

For both cases the scaling law (jl.3|) still holds, and we define backward and stationary 
SS and DSS solutions in the same manner. In particular, a backward SS solution satisfies 
(|1.4|) with a < 0, a stationary SS solution satisfies 

u(x) = X(x)U(X(x)x), X(x) = — , (1.5) 

\x\ 
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with U(x) = u(x), and the profile U(x) of the SS solution for all three cases satisfies Leray's 
equations 

- AU - aU - ax ■ VU + (U ■ V)U + Vp = 0, div£/ = 0, (1.6) 

with a > 0, a < and a = respectively. Note the stationary SS solutions are often called 
minus-one homogeneous solutions in the literature. 

When u(x,t) is either SS or DSS, then so is uq(x). Thus it is natural to assume 

MaOl < rr, O^xgR 3 (1.7) 
| a; | 

for some constant C* > and look for solutions satisfying 

KM)|<^p or |K, tJllis,, < C(C»). (1.8) 

Here by L 9 ' r , 1 < q, r < oo, we denote the Lorentz spaces. In such classes, with sufficiently 
small C* , the unique existence of mild solutions - solutions of the integral equation version 
of (|l.lj) - (|1.2j) via contraction mapping argument, see (12. 5j) - has been obtained by Giga- 
Miyakawa [6] and refined by Cannone-Meyer-Planchon [3111]. As a consequence, if uo(x) is 
SS or DSS satisfying (jl.7p with small C* and u(x, t) is a corresponding solution satisfying 
(11.81) with small C(C), the uniqueness property ensures that u(x,t) is also SS or DSS, 
because is another solution with same bound and same initial data = uq. For 
large C*, mild solutions still make sense but there is no existence theory since perturbative 
methods like the contraction mapping no longer work. 

Alternatively, one may try to extend the concept of weak solutions (which requires uq £ 
L 2 (M 3 )) to more general initial data. One such theory is local-Leray solutions, constructed 
by Lemarie-Rieusset [12] (to be defined below). However, there is no uniqueness theorem 
for them and hence the existence of large SS or DSS solutions was unknown. 

In a surprising recent preprint [7], Jia and Sverak constructed SS solutions for every SS 
uq which is locally Holder continuous. Their main tool is a local Holder estimate of the 
solution near t = 0, assuming minimal control of the initial data in the large (see Theorem 
I3.2p . This estimate enables them to prove a priori estimates of SS solutions, and then show 
their existence by applying the Leray-Schauder degree theorem. Note that this existence 
theorem does not assert uniqueness. In fact, non-uniqueness is conjectured in [7]. 

In this note, as an attempt to understand [7], we consider the existence of discrete 
self-similar solutions for DSS uq satisfying (|1.7p with large C*. 

We now recall the definition of local-Leray solutions, see |12t [7]. Suppose 

u eL| O( ,(l 3 ;l 3 ), II^oIIl 2 := SU P I / \uo\ 2 {x)dx) < oo, divu = 0. (1.9) 

uloc x m 3 KJb^xo) ) 

A vector field u £ Lf oc (M. 3 x [0,oo)) is called a local-Leray solution of ([LT]) - (fL2]) if (i) 

-R 2 



ess 

0<t<i? 2 



sup sup / \u(x,t)\ 2 dx+ sup / / |Vii(x, t)\ 2 dx dt < oo, (1-10) 

<R 2 x <=R 3 JB(x ,R) x gM. s J0 Jb(x ,R) 

and lim [ f \u(x, t)\ 2 dx dt = 0, (1.11) 

\x \^coj JB(x ,R) 
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for any R < oo, (ii) together with some distribution p in they satisfy (jl.ip in in the 
sense of distributions, (hi) lmn_>o+ ll u (">*) ~ u o\\l 2 (k) = f° r an Y compact set C M 3 , and 
(iv) u is suitable, i.e., it satisfies the local energy inequality in the sense of Caffarelli, Kohn, 
and Nirenberg [2]. 

The class of local-Leray solutions contains both Leray-Hopf weak solutions (with uq G 
L 2 (]R 3 )) and mild solutions (with uq in L 3 (M 3 ) or V MO" 1 ), and is strictly larger. It is 
useful for our purpose because it allows initial data of the size |ito(x)| ~ because the 
local energy inequality is valid, and because a priori local energy estimates are available 
(see Lemma [3~T|) . 

We now state our main theorems on the existence of forward discrete self-similar solu- 
tions. We first consider those with DSS factor close to one. We denote (z) = (\z\ 2 + 2) 1 / 2 
for z G W 1 , n G N. 

Theorem 1.1 (Existence of DSS solutions with factor close to one). 

For any < 7 < 1 and C* > 0, there is A* = A* (7, C*) G (1, 2) such that the following 
hold. Suppose uq G C^ c (]R 3 \{0}), ll^oll^y^^ n < C*> divuo = 0, and uq is DSS with 
factor A G (1, A*]. Then there is a local-Leray solution u of (jl.ip with initial data uq that 
is DSS with factor A and, for v(-,t) := u(-,t) — e tA UQ 

\u(x,t)\ < —^=, \v(x,t)\<^- t (1.12) 

in with C = (7,6**). It is also a mild solution in the class f|l . 12 j) i . If furthermore, 
\\ u °\\c 1 ^(b 2 \b 1 ) — ^* f or some < /3 < 1, then 



\v(x,t)\ < 



in R\ with C = (/3, C*) 



u^y 3> °^)- iDMi - t)i <li3) 



Note that A — 1 > has to be small enough. 

A similar result is true for axisymmetric initial data with no swirl that is DSS with 
arbitrary factor. We recall that a vector field u in M 3 is called axisymmetric if in cylindrical 
coordinates r,9,z with (xi,X2,x%) = (r cos 9, r sin#, z) and r = y 'x'f + x\, it is of the form 

u(x) = u r (r, z)e r + u e (r, z)ee + u z (r, z)e z . (1-14) 

The components u r ,u s ,u z do not depend upon 9 and the basis vectors e r ,eQ,e z are 

er= pi,2£, V e e =(-^,^,0), e z = (0,0,1). (1.15) 
\ r r / V r r J 

It is called "no swirl" if u e = 0. This class of vector fields is preserved under (jl.ip . If 
the initial data no G -£f 2 (M 3 ) is axisymmetric with no swirl, global in-time regularity of the 
solution was proved independently by Ukhovskii-Yudovich [20] and Ladyzhenskaya [8]. See 
|llj for a refined proof. The case of general axisymmetric flow with u e 7^ is open. 



Theorem 1.2 (Existence of axisymmetric DSS solutions with no swirl). 

For any 1 < A < 00, < 7 < 1, and C* > 0, suppose u$ G C^ c (M 3 \{0}), is axisymmetric 
with no swirl, DSS with factor A, divuo = 0, and \\uq\\ , „ , < C*. Then there is a 
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local-Leray solution u of (jl.ip with initial data uq that is DSS with factor A, axisymmetric 
with no swirl, and satisfies (| 1 . 12[) in with C = C(A,7, C*). It is also a mild solution in 
the class (fL"T2"|) i . 

If furthermore, || u o|| c .i, J 8/g^y B . l \ — C* f or some < (3 < 1, then it satisfies (11 . 13f) in 
withC = C(A,/3,C*). 

If C* is small, the existence is known by [6j OH]. Theorems 11.11 and 11.21 are concerned 
with large C*. 

If one assumes higher regularity of uq, say no € (R 3 \{0}) for some < (3 < 1, the 
same proof of Th 4.1] shows 

Ct \t\ 

This rate is optimal in view of the explicit spatial asymptotes for small SS solutions with 
smooth initial data in [TJ. Eq. ()1.13|) is slightly worse than ()1.16|) by a log factor, but only 
assumes no G C\f c . There is a gap between (|1.12p and (|1.13p especially when one takes 
1 — 7 = f3 <C 1. It is probably because we require pointwise bound of the source term when 
we estimate the Stokes system. One may be able to narrow the gap by considering integral 
bounds of the source term. 

Our approach follows that of [7j, and relies heavily on the a priori estimates of the 
solutions, see Lemmas 13.31 and 13.41 One difference is that, instead of estimating a stationary 
solution of (jl.6p , we need to estimate a time dependent solution of (jl.ip . Another difference 
is the following: [7] first proves a priori estimates and constructs solutions for smooth initial 
data, and then gets solutions for C 7 data by approximation. In contrast, we prove a priori 
estimates and construct solutions for C 7 initial data directly. The reason for this change is 
that, at least for Theorem ll.il we need the explicit dependence of A* on the local C 7 -norm 
of the data. 

To extend these results, one may look for DSS solutions with DSS initial data of the 
form 

uq = u\ + uI ( 1 - 17 ) 

where u\ is SS and large, while u\ is DSS with a large factor and is sufficiently small. When 
A is large, a priori estimates seem unavailable, and one may try to study the linearized flow 
around u , a SS solution with initial data Uq. It turns out to be very challenging. 

Another interesting problem is the non-uniqueness of local-Leray solutions conjectured 
by [7] and Sverak [18]: Considers SS solutions W a with SS initial data auo, a > 0. For 
a small, W a is unique (see Lemma I3.4p . However, when one increases cr, one might get 
bifurcation. If the bifurcation is of saddle-node type, we get two SS solutions with the same 
initial data. If it is a Hopf bifurcation, the new solutions would be time periodic in the 
similarity variables (y = t~ l l 2 x and s = logi) and correspond to DSS solutions. These 
kind of DSS solutions u are different from those in Theorems 11.11 and 11.21 since their initial 
data uq are SS and only the difference v(-,t) = u(-,t) — e tA uo are strictly DSS. One may 
approximate uo by DSS data Uq and take limits e — > 0, and try to show that the strict DSS 
property of the corresponding solutions u £ is somehow preserved in the limit. Of course 
this is purely speculation. 

The existence question of discrete self-similar solutions also occur in two other instances 
for Navier-Stokes equations: (i) For singular backward solutions u(x, t) : M 3 x (— oo, 0) — > M 3 
of (jl.ip . the nonexistence of SS solutions under some minimal integrability assumptions 
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was proved in |14j and [19]. The existence problem for DSS solutions under the same inte- 
grability assumptions is open; (ii) For stationary Navier-Stokes equations, the self-similar 
(minus-one homogeneous) solutions in ]R 3 \{0} are shown by Sverak |17j to be exactly those 
axisymmetric solutions found by Landau [9] [10] . The existence problem for strictly DSS 
solutions is open. 

The rest of this note is structured as follows: In §2 we consider the Stokes system. In §3 
we prove a priori estimates for DSS local-Leray solutions. In §4 we show their uniqueness 
for small initial data. Finally in §5 we prove their existence. 

Notation. We denote (z) = (\z\ 2 + 2) 1 / 2 for z <E W 1 , n G N, and A < B if there is a 
constant C, which may change from line to line, such that A < CB. We denote by D^u all 
k-ih order partial derivatives of u with respect to the variable x. 

2 Stokes system 

Consider the non-stationary Stokes system in M 3 with a force tensor / = (fij) 

d t v - Av + Vp = V ■ f, dwv = 0, v\ t=0 = 0. (2.1) 
Here (V • f)j = Y2k ®kfkj- If / nas sufficient decay, a solution is given by v = <&/, with 

(§f)i(x,t)= f f d Xk S ij (x-y,t-s)f kj (y,s)dyds (2.2) 
Jo Jm 3 



and S = (Sij), the Oseen tensor, is the fundamental solution of the non-stationary Stokes 
system in M 3 (see [15] and (THJ, page 27]) 



S ij (x,t)=T(x,t)S ij + ]-^- f P^\dy, Q i ( x ,t) = M^, (2.3) 

4vr dxidxj J R3 \x - y\ 4vr \ x P 



where T(x,t) = (47rt)~ n ^ 2 e~ x2 ^ 4t is the heat kernel. It is known in Theorem 1] that 

Did*s(t,x)\ < c k ,i(\x\ + Vty 3 ~ e - 2k , (e,k> o), (2.4) 

where D l x indicates ^-th order derivatives with respect to the variable x. 

When the bilinear operator &(u®v) : X x X — > X is well-defined on some Banach space 
X of M 3 -valued fields on we say u(x,t) is a mild solution of (jl.ip and (|1 .2j) if 

u(-,t) = e tA u - $(n® u)(-,t) in X (2.5) 

We start with an integral estimate. 

Lemma 2.1. Let 0<a<5, 0<6<5 and a + b > 3. Then 



4>(x, a, b) 
is well defined for i£i 3 and 



f [ (\x -y\+ t)- a (|y| + Vt)- b dy dt (2.6) 

Jo Jr 3 



<f>(x, a, b) < R~ a + R~ b + R 3 ~ a ~ b [l + (l a=3 + 1 6=3 ) log R] (2.7) 
where R = \x\ + 2. 
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The cases a, b £ {3, 4} are stated in [3 (4.12)]. We will also use 4 < a < 5 and 2 < b < 3. 
Proof. Clearly eft < 1 for R < 8. Consider now R > 8. Estimate the integral in 3 parts: 

(•)< [ I \y\- a - b dydt = CR 3 - a - b , (2.8) 

'0 J\y\>2R JO J\y\>2R 



f I (•)< f I R- a (\y\ + Vt)- b dydt 

JO J\v\<R/2 JO J\v\<R/2 



10 J\y\<R/2 JO J\y\<R/2 

<R~ a \ C I + C I )(-)<R~ a (l + R 3 ~ b (l + l b=3 logR)), 



J\y\<l JO J\<\y\<R/2 



and similarly 



/ / (•) < I f (k -y\ + VT^~ty a R- b dydt 

10 JR/2<\y\<2R J J R/2<\y\<2R 

<R- b f [ (\z\ + Vt)- a dzdt<R~ b (l + R 3 ~ a (l + l a=3 logR)). 

JO J\z\<3R 



/ /- \ 2+m 

Lemma 2.2. Suppose \f(x,t)\ < \ rf^J inM.% for0<m<l. Then 



\x\+Vt 

( V~t 



2+m 

\*f(x,t)\ < (ur^J > V M) G <■ 

Proof. By (12.4f) and change of variables x = y/ix, y = Viy, s = ts, 

l*/0M)l< ft (\x-y\ + Vt^~s)- 4 s m / 2 (\y\+V~sr^ m dyds (2.12) 
Jo Jrz 

= r 1 ' 2 f 1 [ (\x-y\ + VT^~s)-*s m / 2 (\y\ + V~ S )- 2 - m dyds. (2.13) 
Jo Jm. 3 



(2.9) 



(2.10) 



□ 



By s m l 2 < 1 and LemmaED we get |S/(s,t)| < t~ l l 2 (x)- {2+m) , i.e., (I2~m . □ 

We now show Holder estimates in space and time. Fix < 9 < 1. Denote a local 
parabolic Holder estimate for (x,t) E ]R^_: 



[u] e (x,t) := sup 



\u(x, t) — u(x, t)\ 



x.t 



6:=\x-x\ + x /\t-t\<^\. (2.14) 



Lemma 2.3. Suppose \ f(x,t)\ < (\x\ + y/i) 2 in Then $/ is locally Holder continuous 
in x and t with any exponent < 6 < 1 and /or any T € (1, oo) 

[$f]o(x,t) < C T {x)~ 2 , VxGM 3 , Vl<t<T. (2.15) 



Proof. We may assume £ = 1. 

We first show spatial Holder estimate. For h G M 3 with 5 = \h\ < 0.1, we have 

\<f>f(x + h,t)-<S>f(x,t)\<h+I 2 (2.16) 
■=[([ +[ )\D x S(z + h,s)-D x S(z,s)\-\f(x-z,l-s)\dzds (2.17) 

JO y./|z|>2<5 J\z\<26j 



For Ji, by mean value theorem and (J2T 



By Lemma |2.1| 



h< f [ \h\\D 2 x S(z,s)\-\f(x-z,l-s)\dzds (2.18) 

JO J\z\>2S 

< f 1 [ S e (\z\+^)~ 4 ~ 8 (\x- z\ + VT^)~ 2 dzds (2.19) 

Jo J\z\>2S 

-2 



For I 2 , we have |^ + h\ < 35. If |x| < 45, splitting 0<£<lto0<£<i and ~ < £ < 1 



/i < (2.20) 

»i| < 3(5. If |z| 

and using (|2.4p . we have 

h< [ 1/2 [ {\z\+yflyUzd S + ^ [ (\z\+VT^7sr 2 dzds 
Jo J\z\<35 Jl/2J\z\<7S 



(2.21) 



Using 

L 1 Ik/ 121 + ^ ds 5 { /.og(i/i), < ^4, (2 ' 22) 

(which can be shown by splitting (0, 1) = (0, 5 2 ) U [5 2 , 1)), we get 

h < 5 + <5 3 log(l/(5) <5. (2.23) 
If 45 < \x\, we have (using (|2.22p again) 



h< C f (\z\+^S)- 4 (\x\ + VT^r 2 dyd S 

Jo J\z\<38 



'\z\<36 

< [ 1 I {\z\ + ^)- 4 (x)- 2 dzds+ [ I {\x\ 2 + 1 - s)~ l dzds (2.25) 

JO J\z\<3& Jl/2J\z\<36 

-2 , X 3 



(2.24) 



'g^-^<S(x)- 2 , (2.26) 



which is also bounded by the right side of (|2.20p (and much less if <5 <C 1). 

We next show temporal Holder estimate. Take r = 5 2 with < 5 < 0.1. (For r < we 
can reverse £ and £ + r). We have 

|$/(x, 1 + r) - $/(x, 1)| < /i + / 2 + J 3 (2.27) 
:= y \DS(z, 1 + t-s)- DS(z, 1 - s)\ ■ \f(x -z,s)\ dz ds (2.28) 
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+ J J J \DS(z, 1 + r - s)| • - z,s)\dzds (2.29) 



l-T 



+ / /I - D5(z, 1 - a) | ■ |/(x - z, s)| dz ds. (2.30) 



For 7i, by mean value theorem and ([2 

•l-r 



/i < jf ^Tdzl +VT^7)- 6 (|x-«| + v^)~ 2 ^ds (2.31) 
~Jo ' f 6e (\ z \ + ^^s)~ 4 ~ 6 (\x-z\ + ^)~ 2 dzds. (2.32) 

By Lemma 12.11 

h < S d (x)- 2 . (2.33) 
The two terms I<i and I3 are similar and it suffices to estimate I3: By (|2.4p . 

/ 3 < ^ y (kl + v 7 !^) -4 (|x - z| + l)" 2 dz ds. (2.34) 
Integrating in s first, we get 

T z< I T^—\ z \~ 2 {\x-z\ + l)- 2 dz. (2.35) 

J R 3 \Z\ 2 + T 

If |x| < 2, then |x — z| + 1 ~ (z) and 

h< J ^\z\~ 2 {z)- 2 dz<5 e . (2.36) 

If Ixl > 2, 



J 3 < / ^ \z\~ 2 \x\~ 2 dz + t\z\~ 4 \x - z\~ 2 dz (2.37) 

■/|2|<M/2 \Z\ + T y|«|>|a;|/2 

= C5\x\~ 2 + C'(5 2 |xr 3 £ ^kr 2 - (2-38) 
The equality here is obtained by rescaling. □ 

Finally we give a Liouville lemma. 

Lemma 2.4. Ifv(x,t) : R\ ->■ M 3 satisfies \v(x,t)\ < Ct~ 1 / 2 (x/y/t)~ 1 ~~ l for some < 7 < 
1 and 

<9 4 v - + Vp = 0, div « = 0, v\ t=0 = 0, (2.39) 
/or some distribution p, then v = 0. 

It is similar to Lemma 4.1 (i)], with exactly the same proof. 



3 A priori estimates for DSS solutions 

We first recall a couple estimates for local-Leray solutions from [7]. 

Lemma 3.1 ([7j Lemma 3.1). There are constants < C\ < 1 < C2 such that the following 
holds. Suppose divuo = 0, A = sup XogK 3 Jb r ( Xo ) \ u oi x )\ 2 dx < 00 for some R > and u is 
a local-Leray solution with initial data uq. Then for A = C\ min {A~ 2 R 2 ,l), 



esssup sup / \u(x,t)\ 2 dx+ sup / / \Vu(x, t)\ 2 dxdt < C2A. (3.1) 

0<t<XR 2 x £R 3 J B R (x ) x m 3 Jo Jb r {x ) 

Because (jl.lip holds for u, for some p(t) = p Xo ,ii(t), 

f XR 2 r 

sup / / \p(x,t) -p{t)fl 2 dxdt < C 2 A Z I 2 R 1 / 2 . (3.2) 



Theorem 3.2 Q7J Th 3.2). Suppose div^o = 0, ||uo||?2 < A < 00, and \\uo\\ci(B2) — 

uloc 



M < 00 for some 7 £ (0, 1) . T/ien i/iere exists T = T(A, 7, M) > such that any local-Leray 
solution u with initial data uq satisfies u G Cpar(-Bi/4 x [0, T]) and 

ll^ll^^xp.TD^^^M). (3.3) 

We now consider DSS solutions with factor close to one. 

Lemma 3.3 (A priori estimates for DSS solutions with factor close to one). 

(i) For any < 7 < 1 and C* > ; there is A* = A* (7, C*) € (1,2) such that the 
following hold. Suppose 1 < A < A* and u is a forward X-DSS local-Leray solution of (jl.ip 
with X-DSS initial data uq G C^ c (]R 3 \{0}) satisfying divuo = and \\ u o\\ c ^^§^ Bi ^ < C*- 
Then v(x,t) := u(x,t) — (e tA uo)(x) satisfies, for some C = C(j,C m ), 

\u{x,t)\ < ° \v(x,t)\ < ^=(^=\ , (x,t)£R%. (3.4) 



x 



Moreover, u is a mild solution of (jl.l|) and (|1.2|) . 

(ii) Lf furthermore \\ u o\\qi,p(b 2 \ b s < C* /or some < /3 < 1, then 

m /or some C = C(/3, C*). 
Proof, (i) We will first show a weaker estimate 



(3.6) 



We first consider the region below the paraboloid, 

t<\x\ 2 /R\, (i,t)g^, (3.7) 
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for some Ri = i?i(7, C*) > sufficiently large to be decided later. By Theorem 13.21 there 
exists T\ = Ti(7, C*) > such that for any xq G M 3 with 1 < |xq| < A 



W^vWJ^Bygixo) x [0,Ti]) < C( 7 ,a), (3.8) 
where we have used that e* A iio satisfies the same Holder estimate. Since v (x, 0) = 0, we get 

\v(x,t)\ < C7( 7 ,a)^/ 2 , I < |^| < A + 0<£<7i. (3.9) 

Since \ k G Z, is another local-Leray solution with same initial data uo, the above 
estimate remains valid with v replaced by z/ A Scaling back, 

H X ,t)\< C ^ /2 , (3.10) 

in cylinder C k = {x : §A fc < \x\ < X k+1 } x [0, \ 2k Ti) for every fceZ. Since § < A < 2, the 

union U^g^Cfc contains a set of the form in (|3.7p with iii = 2T X 

For the complement, by rescaling it suffices to prove an upper bound in the region 

t>\x\ 2 /Rj, l<t<\ 2 . (3.11) 

This region satisfies |x| < XR±. Let T = 4. (We will take T = 4A 2 for the proof of Lemma 
El) By LemmaEH u G {LfL 2 x n L 2 t Hl) C lJ° /3 in B XRl x (0, T) with the norm bounded 
by C(7, C*). Together with the decay for e tA uo and estimate (|3.10|) for v in the region (13. 7p . 
we get ll'"ll jL io/3 (R 3 x(0 ^ )) < C(7,C*). Since p = (-A)~ l didjUiUj, we get 

iHi<e (K 3 X(0 ,T))^^^)- ( 3 - 12 ) 

Denote shrinking annuli A fc = {x G R 3 : 2Ai?i + k < \x\ < 2XR 1 + 20 — k} for k = 1, 2, . . .. 
Note that -u G L°°(Ai x [± T]) since the region is contained in the region (|3.7p . By regularity 
theory for D x u G L°°(A 2 x [1,T]) for £ < 3. By -Ap = {diU^djU,) and (I3TT21 . we 

get i/p G L°°(A 3 x [1,T]) for i < 2. By (JTIJ, L>V G L°°(^ 3 x [1,T]) for l<\. 

Choose a smooth cutoff function £(x) > with C,{x) = 1 when \x\ < 2\R\ + 5 and 
Q{x) = when \x\ > 2\R\ + 6. Let w(x,t) be a solution supported in x G of 

div«;(x,t) = u(a;,t)-VC(x), IK,*)||^ < IK,*)!!^^) VI < t < T, (3.13) 

by a construction uniform in t, see e.g. [5| §111.3]. In particular divc^ w = dtu ■ V£ and 
\\d t w(-,t)\\ H i < \\d t u(-,t)\\ L 2 [A3) . Let 

fi = £u-u;. (3.14) 
One can check that u is a suitable weak solution of 

d t u - Au + (u ■ V)u + Vp = /, divu = 0, (3.15) 

satisfying 

u(x,t) = if |x|>2Ai?i + 6; /(x,t) = if x g A 4 , (3.16) 

lL t °°Hi(R 3 x(l,T)) < C(l,C*), (3-17) 



10 



and by Lemma 13. II 

ess sup j \u(x,t)\ 2 dx+ [ [ \Vu\ 2 dx dt < C(-y, C*). (3.18) 
0<t<T JR3 Jo y R 3 

By d3H|), there exists a time ti G [1,3/2] such that u(-, t{) G i^QR 3 ) with ||£t(-,ii) < 
(7(7,(7*). Thus there is T 2 = T 2 (7, (7*) > and a strong solution u(x, t) : M 3 x [ii, ti + T 2 ) — > 
M 3 of (|3.15p with initial condition u{x,t\) = u(x,t\). We may assume T 2 < 1. 

By weak-strong uniqueness, we have u(x,t) = u(x,t) for (x,t) G R 3 x [t\,ti + T 2 ). 

If we take t 2 = h + T 2 /4 < 2 and A 2 = 1 + T 2 /4, we have [t 2 ,t 2 X 2 ] m (h,h + T 2 ) 
whenever ii G [1,3/2] and 1 < A < A*. All spatial derivatives of u are a priori bounded in 
M 3 x [t 2 ,t 2 A 2 ]. Because u is discrete self-similar and agrees with u in the region (|3.1ip . we 
get \u(x, t)\ < C in the region (I3.11j) . 

We have shown the weaker estimate (|3.6p . We now show (|3.4p . Let / = —u ® u. By 
|(e' A u )(x)| < r 1 /2(i- 1 /2 s >" 1 and (1331) . 

I/M)I<-TZ7> (M)GR£. (3.19) 

2T + E 

By Lemma [2.21 (with m = 0), we get \$>f(x,t)\ < r _1 / 2 (x/v / i) 2 , where $ is defined in 
(|2.2p . Note v = v — &f satisfies the bound (|3.6p and the linear Stokes system in with 
zero initial data and zero source. By Lemma 12.41 v = 0. Thus we have v = $/ and (13.4p . 

(ii) Suppose now no is in and ||V-uo|| g/3 ^^ Tli) ) — ^C* for any xq G IR 3 with 
1 < |^o I < A. The vorticity u) = curlu satisfies 

d t uj - Auj = - curl V • (u <g) u), (3.20) 

and w = w(-,0) satisfies ll^ol| C;9( g(a;0]1/2 ^ < CC*. 

Denote = | — and Qfc = B(xo,rk) x (0, Ti] for k = 0, 1, The previous step 

shows u G Cp 0r (Qo)- 

Let r/(x) = 7/0(2; — xo) where t]q{x) is a fixed smooth cut-off function with t]q(x) = 1 for 
|x| < r 2 and r]o(x) = for |x| > r\. Decompose u) = uj\ + u 2 + where 

= e A ^ s )[curlV- (it ® wj)](.,*)ds, 
Wa (. J t) = e* A (i^(- > 0)), (3 ' 21) 

W3 = U) — UJ\ — Ld 2 . 

Note (dt — A)i(j3 = in Q 2 and W3(a;,0) = for x G B r2 (xo). Thus, with possibly a smaller 
T\, both u> 2 and u>3 are bounded in Cp ar (Q3) by (7(/3,(7*). By singular integral estimates 
for heat equation, we have 

IMIz/HRSxlo.T!]) < C||«®«»;IU«(RSx[0,2\]) < C(P,C*,q) (3.22) 

for any 1 < q < 00. We take g = j-^. Thus uj = Ylf=i w « ^ L q (Qz). By elliptic estimate, 

||V«|| W (q 4 ) < C\\cmlu\\ Lq (Q 3 )+C\\divu\\ Lq (Q 3 ) + C\\u\\ Lq (Q 3) < <7(/3,C*). (3.23) 

We now do a similar decomposition of u with 770(2:) = 1 for \x\ < r@ and 7/0(2;) = for 
\x\ > 7-5. Again uj 2 and W3 are in CpariQl)- Rewrite 

V ■ (u®uri) = u- V(u? ? ) G L g (R 3 x [0, Ti]). (3.24) 
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By heat potential estimate 

[wi] c ^ r(R 3 x[0iTl]) < C\\u- V(n?7)|| L9(IR 3 x[0iTl]) < C(/3,C*). (3.25) 

Thus uj = Yli=i w « e Cpar{Qi)- By elliptic estimate, 

\\ Vu \\ L ™^cnW^)) - C\\cuTlu\\ + C\\d\vu\\+C\\u\\ < C(P,C.) (3.26) 
where the middle norms are Cp ar (Q7)-norms. In particular we have shown 

\Vu(x,t)\ < C(fi,C*), (1 < |x| < A, 0<t<Ti). (3.27) 
By the same scaling argument for (13. 10f) . we get 

\Vu(x,t)\< , C ^'°i 2 (3-28) 

(Vt + \x\Y 

in the sub-paraboloid region ([3 
We now rewrite 



with 



Vi(x,t)= / Sij(x - y,t - s)gj{y,s)dyds (3.29) 



g = -u-Vu, \g(y,s)\ < ,^?*L (3.30) 
Thus for i = 0, 1, by (12. 4|) and change of variables x = y/ix, y = yty, s = ts, 

\D e v(x,t)\< [ [ (\x-y\ + Vt^)- 3 - e {V^+\y\)~ 3 dyds (3.31) 

JO JM 3 

= t -(l+/)/2 f 1 f ( |^_ y | +v ^ ) -3-^ ( ^ +M) -3 dyds _ (332) 

By Lemma [27X1 we get ([33]). □ 

We next consider axisymmetric DSS flow with no swirl. 

Lemma 3.4 (A priori estimates for axisymmetric DSS flow with no swirl). 

(i) For any 1 < A < oo, < 7 < 1, and C* > 0, suppose u is a forward X-DSS local- 
Leray solution of (jl.ip with X-DSS initial data uq € C^ C (1R 3 \{0}) that is axisymmetric with 



no swirl, DSS with factor \, divun = 0, and I kin I Lw"5 \ D \ < C*- Then u satisfies (137 

( ^ T (.- fc >A\-t>lJ 

with constant C = C(A,7, C*). Moreover, u is a mild solution of (jl.ip and (|1.2p . 

(mJ If furthermore \\uo\\ cl ,g * < C* /or some < /3 < 1, i/ien (|3.5p /ioW ttratft 
constant C = C(A, /3, C*). 

Proof. We use the same proof of Lemma [3731 (i) until we get a time ii € [1,3/2], a T2 = 
^2(7, C*) E (0, 1) and that u agrees with a strong solution in [t\,t\ +T2). For £2 = ti +T2/4, 
all spatial derivatives of u(x,t2) are a priori bounded. Since tt has compact spatial support, 
we get 

||«(-,t2)||fl3(R3) <C(7,a). (3.33) 

By [HI Th 1], 

\\u(;t)\\m(R 3 ) < C(7,a), ||n(.,t)||^ 2(R 3) < C( 7 ,C*,t), (3.34) 

for t2 < t < T = AX 2 , which contains the interval [£2^2 A 2 ]. Since ii is A-DSS and u agrees 
with u in the region t > \x\ 2 /R\ and 1 < t < T, we have shown the boundedness of u in 
the same region. The rest of the proof of Lemma 13.31 then goes through. □ 
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4 Uniqueness of DSS solutions with small data 

When the initial data no is small in L 3 '°°(IR 3 ), the existence theorem of [6J [3j d] says that 
there is a unique mild solution u m -,^ (x, t) in the class 

ll M mild(*)||BC„([0,oo);L3,oo( R 3)) < C||uo||i3,<x> ( R 3) . (4.1) 

Above BC W means bounded weak-star continuous L 3,00 -valued functions of time. It is also 
known that n m iid is a local-Leray solution, see [12]. However, for our application later, we 
will need uniqueness in a larger class of solutions. 

Lemma 4.1. Let uq satisfy the assumptions of Lemma \3.3\ or Lemma \3.4\ with C* suffi- 
ciently small. Then any X-DSS local-Leray solution u(t) of (jl.ip with initial data uq must 
agree with the mild solution vtimiid 

constructed by JH [5j \J$. 

In particular, u(t) is allowed to be large in BC w ([0, oo); L 3 '°°(]R 3 )). Such a statement 
that "large equals small" is not known for general solutions with small _L 3 '°° data. The 
lemma is only for DSS solutions and relies on the estimates of Lemmas 13.31 and [37 



Proof. Denote W(x,t) = u m nd(x,t) and v(x,t) = u(x,t) — u m iid(x, t). They are both DSS 
with factor A and satisfy (|3.4|) . Thus 

|W ' (M)I < (RT7I)' |l,(l ' t)l < (kf+W' (4 ' 2) 

in Note that v satisfies 

d t v - Av + (W + v) ■ Vv + v ■ VW + Vp = 0, divv = 0. (4.3) 

We have -Ap = Yjij didj((W+v)iVj+ViWj) and hence \\p(-,t)\\ L ^ R 3^ < || (x)~ 3 || L9 ( R 3) < oo 
for 1 < t < A 2 and 1 < q < oo. 

Let Cr( x ) = C( x /R) where ((x) > is a smooth cut off function with ((x) = 1 for 
\x\ < 1 and £(x) = for \x\ > 2. Multiplying (|4.3|) by vQr and integrating by parts over 
M 3 x [1,A 2 ], we get 

[ ^ X ^ 2 ( R (x)dx\ + [ X [ {\Vv\ 2 C R -v®W: (Vv)Cr) dxdt = L r (4.4) 
Jr3 2 J t= i J i Jm.3 

where 

T R -= I I \\v\ 2 ACr+(}^-{W + v) + {W -v+p)v)-VC R dxdt. (4.5) 
Ji Jm.3 2 2 

By Lemma 13.31 and p £ LfL%, I R converges to as R — > oo. The first term in (I4.4|) also 
converges and it converges to 



\v(x,t)\ 2 



dx 



x 2 

- A ~ 1 
t=i = 



\v(x, l)\ 2 dx > 0. (4.6) 



We have shown 

■A 2 



[ [ {\Vv\ 2 Cr-v®W : (Vu)Cfl) dxdt<o{\) (4.7) 
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as R — > oo. Since 

v ® W : (Vv)6l dxdi < i ^ \Vv\ 2 (r + C J J \vW\ 2 Cr (4.8) 

and the last term converges as R — > oo, we get JJ \Vv \ 2 (r Jf |Vt>| 2 < oo as R — > oo and 
by Hardy inequality 

\Vv\ 2 < JJ v ®W :Vv < cJJJ i^fY [JJ \Vv\ 2 Y < C* J J \Vv\ 2 . (4.9) 

Thus when C* is sufficiently small we get /j* / R3 |V-y| 2 ctecft = 0. Thus u = 0. □ 

5 Existence of large DSS solutions 

In this subsection we prove Theorems 11.11 and 11.21 Let 

U(x,t) = (e tA u )(x). (5.1) 
By the assumption on tio, U is A-DSS and 

Introduce a parameter a £ [0, 1]. We look for a solution u(x,t) of (|1.1[) of the form 

u(x, t) = aU(x, t) + v(x, t), u(x, 0) = uuq{x). (5-3) 
The difference v satisfies the nonhomogeneous Stokes system (|2.1j) with 

f = -(aU + v)®(aU + v). (5.4) 
We expect (which is true at least for small a) that v is A-DSS and 

f{x,t) = \ 3 f(Xx,\ 2 t), V(z,t)€R|. (5.5) 
In view of Lemmas 13.31 and 13.41 we also expect 



The decay rate of v can be improved but it is not needed. 

We now set up the framework for the application of the Leray-Schauder theorem. Let 

Q = R 3 x[l,A 2 ], (5.7) 

and the Banach space X = X(X): 

{, divv = 0, \\v\\x < oo, 
vGC(Q;R 3 ): " , (5.8) 

v(x, 1) = Au(Ax, A ), Vx G IR J J 
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where 



\\v\\x '■= sup (x) ' y \v(x,t)\. 
(x,t)eQ 

For each v G X, we define its DSS extension by 



(5.9) 



Ev(x,t) = X k v(X k x,X 2k t), for (x,t) G H.| 



(5.10) 



where /c is the unique integer so that 1 < X 2k t < A 2 . 
We now define an operator K : X x [0, 1] — > X by 



K(v, a) := -$[(ctC7 + (af7 + Ev)]\q, \/v G X, Vex G [0, 1]. 



(5.11) 



Above $ is defined by ([23]) . 

Note that for v G X with < M and < a < 1, the force / = -(erf/ + EJu) 8> 

(<j£7 + Ev) satisfies ([53]) and and $(/) defined by (Jl^D satisfies dHJ) and is A-DSS. 

By Lemma [2.21 its restriction to Q, K(v,a) = — $(/)|q, is inside X and ||^(v, o")||x < 
C(C* + M) 2 . Thus iiT indeed maps bounded sets in X x [0, 1] into bounded sets in X. 

Furthermore, K is compact because its main term §{aU ctU)\q is one dimensional 
while the other terms of K have extra decay by Lemma 12.21 and are Holder continuous in x 
and t by Lemma 12.31 

We have now a fixed point problem 



that satisfies the following: 

1. K{v, a) : X X [0, 1] — )■ X is compact by the previous discussion, 

2. it is uniquely solvable in X for small a by O El H] and Lemma |4.1| thus the Leray- 
Schauder degree is nonzero, and 

3. we have a priori estimate in X for solutions v of (|5.12p uniformly for all a G [0, 1] by 
Lemma 13.31 or 13.41 

By Leray-Schauder degree theorem (see e.g. |13j), there is a solution v G X of (|5.12p 
with (7 = 1. It follows that Ev satisfies (|2.ip with / = — (U + Ev) <S> (U + -E'v), and hence 
■u = U + -Eu is a A-DSS local-Leray solution of ([Lip with initial data no- □ 
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